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We study coherence in a cold atom single photon switch where the gate photon is stored in a 
Rydberg spinwave. With a combined field theoretical and quantum jump approach and by employing 
a simple model description we investigate systematically how the coherence of the Rydberg spinwave 
is affected by scattering of incoming photons. With large-scale numerical calculations we show how 
coherence becomes increasingly protected with growing interatomic interaction strength. For the 
strongly interacting limit we derive analytical expressions for the spinwave fidelity as a function of 
the optical depth and bandwidth of the incoming photon. 
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I. INTRODUCTION 

Cold gases of Rydberg atoms are currently receiv¬ 
ing a growing attention in the communities of quantum 
optics m, quantum information (T, and many-body 
physics (SHE]. This is rooted in the fact that they of¬ 
fer strong and long-ranged interactions and at the same 
time grant long coherent lifetimes. Currently, consider¬ 
able efforts are devoted to developing all-optical quantum 
information protocols pH ng with the Rydberg-atom- 
mediated interaction between individual photons [H6UT8] . 
Fundamentally important optical devices that operate 
on the single photon level, such as phase shifters [19], 
switches m and transistors [2T| [22] , have been demon¬ 
strated experimentally in Rydberg gases. 

Single photon switchs might form a central build¬ 
ing block of an all-optical quantum information proces¬ 
sor [23H25] . The prime function of such switches is to 
control the transmission of an incoming photon through 
a single gate photon. One promising way to realize this 
is to store the gate photon in form of a gate (Rydberg) 
atom immersed in an atomic gas which is in a delocal¬ 
ized spinwave state [26ff28] . The gate atom then prevents 
transmission of incident photons through the gas, while 
ideally the coherence of the Rydberg spinwave state is 
preserved j29ff32] . The latter property would permit the 
subsequent coherent conversion of the Rydberg spinwave 
into a photon which would pave the way for gating the 
switch with superposition states that can also be subse¬ 
quently retrieved. Currently, there is only a basic under¬ 
standing of how the coherence of the Rydberg spinwave 
might be affected by the scattering of incoming photons 
and no systematic study of this important question ex¬ 
ists. 

In this work we address this outstanding issue within 
a simple model system. We study the propagation of a 
single photon under conditions of electromagnetically in¬ 
duced transparency (EIT) in a cold atomic gas in which 
a gate photon is stored as a Rydberg spinwave. An in¬ 
cident photon subsequently experiences a Rydberg me¬ 
diated van der Waals (vdW) interaction with this stored 
gate atom which lifts the EIT condition and renders the 
atomic medium opaque. In this case the incident pho¬ 



FIG. 1. (a) EIT level scheme. The groundstate 11), excited 
state 12) (decay rate 7 ) and Rydberg state |3) are resonantly 
coupled by a single photon field £(z,t) (with collective cou¬ 
pling strength g) and a classical field of Rabi frequency O. 
Initially a gate photon is stored as a spinwave in the Ryd¬ 
berg state 14) (indicated by the green circle). (b,c) Polariza¬ 
tion profiles Pj(z,t) for a spinwave consisting of two possible 
gate atom positions Zj (j = 1, 2) and their dependence on 
the blockade radius Rb and the system length L. (b) For 
L > Rb and \Z2 — Z±\ > 2Rb the polarization profiles asso¬ 
ciated with the two gate atom positions are distinguishable, 
(c) When L < Rb the polarization profile is independent of 
the gate atom position which leads to enhanced coherence of 
the stored spinwave. 


ton is scattered incoherently off the Rydberg spinwave. 
We study the photon propagation and explore the de¬ 
pendence of Rydberg spinwave coherence on the interac¬ 
tion strength (parameterized by the blockade radius Rb), 
the system length L and bandwidth of the incident pho¬ 
ton pulse. Our findings confirm that strong absorption, 
i.e. high gain, can be achieved already for large systems 
(L > Rb) while coherence of the spinwave is preserved 
only for sufficiently strong interactions, i.e. L < Rb. In¬ 
tuitively, this can be understood by regarding the scat¬ 
tering of the incoming photon as a measurement of the 
position of the gate atom. When L < Rb this measure¬ 
ment is not able to resolve the position of the excitation 
and hence coherence of the Rydberg spinwave is main¬ 
tained. Our study goes beyond this simple consideration 
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by taking into account propagation effects, a realistic in¬ 
teraction potential and a finite photon band width. The 
results can therefore be considered as upper bounds for 
the fidelity with which a Rydberg spinwave can be pre¬ 
served and re-converted into a photon in an experimental 
realization of a coherent cold atom photon switch. 

The paper is organized as follows. In section II, we 
introduce a one-dimensional model system to study the 
propagation dynamics of single source photons in the 
atomic gas prepared in a Rydberg spinwave state. In 
Sec. Ill, the model system is solved numerically with re¬ 
alistic parameters. We identify the working regime for 
a single photon switch where the source photon is scat¬ 
tered completely. In Sec. IV, we numerically study the 
fidelity between the initial spinwave state and the final 
state after the source photon is scattered. Our calcula¬ 
tion shows that the coherence of the spinwave is preserved 
when L > R\) while the final state becomes a mixed state 
when L < R]j. In Sec. V, We provide analytical results for 
a coherent single photon switch (L > Rb)- We reveal that 
the transmission and switch fidelity depend nontrvially 
on the optical depth and bandwidth of the source photon 
field. We summarize in Sec. VI. 

II. THE MODEL SYSTEM 

Our model system is a one-dimensional, homogeneous 
gas consist ingof N atoms, whose electronic levels are 
given in Fig.jRi. The photon field £(z,t) and the EIT 
control laser (Rabi frequency Q) resonantly couple the 
groundstate |1) with the excited state |2) and |2) with 
the Rydberg state |3). Following Ref. [34], we use polar¬ 
ization operators P(z, t ) and S(z , t ) to describe the slowly 
varying and continuum coherence of the atomic medium 
11)(21 and |1)(3|, respectively. All the operators 0(z,t) = 
{£(z, t ), P(z, t ), S(z, £)} are bosons and satisfy the equal 
time commutation relation, [0(z, £), (T(V, t)] m S(z-z'). 
Initially, the atoms are prepared in a delocalized spinwave 
state with a single gate atom in the Rydberg state |4), 

where k is the wavenumber of the spinwave and | Zi) = 
|li... 4^... Ijv) abbreviates many-body basis with the 
gated atom located at position Z{ and the rest in the 
groundstate. The Rydberg spinwave state is created rou¬ 
tinely in experiments [2014221 35, 36 . When interacting 
with the incoming single photon, the general many-body 
state of this one-dimensional system is expanded as |2j 


£ = 1 at any moment. We have defined 0(z,t) = 
(0(z,t)), i.e. the expectation value of the operator 
0(z,t). Specifically one finds that E(z,t) is the prob¬ 
ability amplitude in the one photon state, P{z , t) and 
S(zR) are the amplitude of one atom in the |2) and |3) 
state, respectively. 

In order to develop a first intuition for the physics 
at work we first consider a spinwave that is delocalized 
merely over two atoms embedded in the atomic cloud 
(see Fig. [l]3,c). We assume furthermore that the inter¬ 
action between atoms in state |3) and the gate atom is 
infinite for distances smaller than the so-called blockade 
radius Rb and zero otherwise. Outside the blockade re¬ 
gion, the photon propagates (along the direction) as 
a dark-state polariton by virtue of EIT [3T . Inside the 
blockade region the medium behaves like an ensemble of 
two-level system. Here the incoming photon is building 
up a non-zero polarization P(z, £), whose modulus square 
is the probability density distribution for finding an atom 
in the decaying state |2) according to Eq. 0 133 • Even¬ 
tually, this leads to the loss of the incoming photon and 
makes the medium opaque. In order to understand how 
such photon scattering affects the coherence of the prop¬ 
erties of the spinwave one needs to analyze the shape of 
the polarization profile. As shown in Fig. [TJd this in gen¬ 
eral depends on the position of the gate atom when the 
system length is larger than the blockade radius Rb < L. 
Here, since L > 4Pb and \Z 2 — Z\\ > 2Pb> it is possible to 
distinguish the profiles Pj{z , t) which are associated with 
the two possible positions of the gate atom. Conversely, 
the polarization Pj(z, t) becomes independent of the gate 
atom position when L < Rb (see Fig. lip). In this case — 
as discussed in detail later — the coherence of the spin- 
wave will be preserved as one can not distinguish gate 
atoms from the scattered photon. 

Let us now consider the actual photon propagation to¬ 
gether with a realistic interaction potential. The dynam¬ 
ics of the system follows the master equation [34] [38] 

p{t) = +7 f dzP(z,t)p{t)P\z,t), ( 2 ) 

Jo 

where the first term on the right-hand side (RHS) is the 
evolution of p(t) under the effective Hamiltonian H e = 
H p + H ap + H a , and the spontaneous decay (with rate 
7 ) from the state |2) is governed by the second term. In 
the effective Hamiltonian, the photon propagation in the 
medium is governed by the Hamiltonian 

H p = —c / dz£\z,t)d z £(z,t). 


I’M*)) = 


+ 


€ + / dzE(z,t)£\z,t) + j dzP(z,t)P\z,t ) 
IdzS(z,t)S\z,i) |M0)>, (1) 


where £ is probability amplitude of the initial spinwave 
state. In the weak field approximation, we will assume 


with the vacuum light speed c. The atom-photon cou¬ 
pling is described by 


#ap = 


-J 


dz 


y P f (z, t)P(z, t) + g£(z, t)P ] (z, t) 


+ n&(z,t)P(z,t)+h. c. 
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where g = y/~Ng s with g s being the single atom-photon 
coupling strength. The vdW interaction between an 
atom in the state |3) and the gate atom at position Zi is 

^ = dzS\z,t)Vi{z)S{z,t). 

The interaction potential depends on the gate atom po¬ 
sition, 

V i (z) = V i (z)\Z l )(Z i \, 

where Vi(z) = C§j(Zi — z ) 6 gives the vdW interaction 
with Cq being the dispersion coefficient. 

For the case of a single incoming photon which we con¬ 
sider here the solution of the master equation © is [38] 

L t 

p(t) = e~ iff ° t p i e ^*+7 J j dzdt'P(z,t')p i P\z,t'), 

0 0 

(3) 

where p\ = |4/at(0))(4/at(0)| and P(z,t) = 

e' l ^ et P(z 1 0) e ~i^ et . The first term on the RHS de¬ 
scribes the unhindered photon propagation through the 
medium, while the second term accounts for the photon 
scattering, i.e. photon-loss from the medium. 


III. TRANSMISSION OF THE SOURCE 
PHOTON 

To calculate © we first treat the dynamics under the 
effective Hamiltonian in the Heisenberg picture. To this 
end we obtain the equation of motion for the expectation 
values 0(z,t) from the corresponding operator Heisen¬ 
berg equation |2T9 . Note, that due to the linearity of 
the equations we can moreover calculate the expectation 
value for each component | Zj) of the Rydberg spinwave, 
i.e. each of the possible positions of the gate atom, sep¬ 
arately. This yields the set of equations 

d t £j(z,t) = -cd z £j(z,t) +igPj(z,t ), (4a) 

d t Pj(z,t) = ~ Pj(z,t ) +ig£j(z,t) + iQSj(z,t), (4b) 

d t Sj(z,t ) = — iVj(z)Sj(z,t ) + iQPj(z,t). (4c) 

where the index j labels the respective spinwave com¬ 
ponent. Alternatively, these equations can be obtained 
from a Heisenberg-Langevin approach HZ]. We solve the 
coupled equations © through a Fourier transform yield¬ 
ing the formal solution for the polarization 

/ +oo 

dcjXj(z)S 0 (uj)e- iujT+i c So dz'xj(z’)' ( 5 ) 

-OO 

Here we have abbreviated T = t — z/c and introduced 
the electric susceptibility 

,. (2| = „_ “~ V A Z ) _ 

XA 1 "n 2 -[w-v;( s )](u + >7/2)' 



FIG. 2. (a) Photon transmission T as a function of the pulse 
duration r for g — IOOO 7 (squares) and g — 5 OO 7 (circles). 
The medium becomes transparent when 77 < 1, i.e. the band 
width Acj of the pulse is large, (b) Photon transmission as a 
function of the coupling constant g for Rb = L (squares) and 
Rb — L/2 (circles). The solid curve is the analytical result 
obtained from Eq. 0. The dashed curve is plotted as a guide 
to the eye. Note, that T « 0 when g — IOOO 7 for both Rb — L 
and Rb = L/2. The data is calculated for rubidium atoms 
with the parameters, L — 20 /mn, H — 27 , and 7 « 2tt x 5.7 
MHz. The blockade radius can be changed through selecting 
different Rydberg states. For example, Rb — L — 20 gm when 
|3) = 11275) and |4) = 11305), where C 6 «4.2xl0 6 GHz/un 6 . 


From Xj( z ) one can actually extract the blockade radius 
as the critical distance at which the vdW interaction and 
the control laser are equally strong. This yields Rb = 
|7C 6 /2tt 2 |V6 m 

The polarization © depends on the Fourier transform 
£o(uj) of the photon field at position z — 0. To be specific 
we take the photon pulse to be a Gaussian at t = 0 which 
is normalized in space, 


£(z, t = 0 ) 


1 

/ 2 2M/4 eX P 

( 7 rc^r 2 ) 1 / 4 


(z-zq) 2 ~ 
2 c 2 t 2 


Here r is the temporal duration of the pulse and zo is 
the initial central position (zq —cr). The band width 
of the pulse is then given by Acj = 1 /r. Note, that it is 
generally not possible to evaluate the formal solution ([5| 
analytically. Moreover, numerical calculations are chal¬ 
lenging since the involved time and length scales span 
several orders of magnitude m ■ 

Let us now calculate the photon transmission as a 
function of the pulse duration r, which to our knowl¬ 
edge has not been examined previously. We de¬ 
fine the transmission of the photon pulse as T = 
J 0 °° dt\£(L, t)| 2 / J 0 °° <R|£(0, t)\ 2 . In Fig. [ 2 r, we show T as 
a function of the pulse width for two values of the atom- 
photon coupling strength g. For fixed pulse length r, we 
find that stronger couplings generally are accompanied 
by a lower transmission. Furthermore, we observe that 
the transmission increases with decreasing pulse duration 
r. This is due to the fact that the pulse contains increas¬ 
ingly more weight on frequency components, which are 
outside the absorption window of the medium. For the 
purpose of complete photon scattering, one thus has to 
utilize narrow frequency band pulses. 
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Next, we briefly discuss the dependence of the trans¬ 
mission T on the strength of the atom-photon coupling g. 
Fig. [ 2)3 shows data for two choices of the blockade radius, 
R]j = L and R\> = L/2. As expected, T decreases with 
increasing g. However, for the system parameters chosen 
here, there is virtually no dependence of T on the value 
of the blockade radius when g = IOOO 7 , where T « 0. 
These findings indicate that one reaches the strong scat¬ 
tering regime when g 7 and Acj <C 7 . This is the 
working regime for the single photon switch where the 
medium becomes opaque for the incident photon. 


IV. FIDELITY BETWEEN THE INITIAL AND 
FINAL STATE 


Focusing on this regime, our next task is to investi¬ 
gate how the photon scattering influences the Rydberg 
spinwave. We quantify the difference between the ini¬ 
tial Rydberg spinwave p\ and the final state pf by the 
fidelity [41: 


F = 


Tr|v7\7f| 


As the initial spinwave is a pure state, this sim¬ 
plifies to F = (1/7V 2 ) Y^j k A jk, where A jk = 

7 / 0 L Jo°° dzdrP^(z,T)Pj(z,r). This shows that a high 
fidelity can be obtained only if the polarization pro¬ 
files Pj(z , r) for each spinwave component are essentially 
equal: Only when Aj k ~ 1 and thus Yhjk Aj k ~ N 2 the 
fidelity is close to one. This is the formal version of the 
intuitive statement that we made earlier in conjunction 
with the discussion of Fig. [l]3,c. 

For completeness we provide a numerical example for 
which we choose R b = L/2 and select only two compo¬ 
nents of the spinwave, where the gate atom is located 
at either Zi = 0 or Z{ = L. The resulting polariza¬ 
tion profile \P(zR)\ 2 is shown in Fig. |3^, b. For Zi = 0, 
non-vanishing polarization is built up within the block¬ 
ade region as long as the photon is inside the medium 
(Fig#). Integrating over time we obtain the intensity 
I p (z) = J 0 °° dr |P(z,t)| 2 which clearly shows a decay to 
zero within a blockade distance Rb (see Fig. 3]:). In con¬ 
trast, for Zi = L, appreciable polarization is built up also 
outside Rb and the profile is peaked at approximately 
z = L — R\) (Fig. [ 3 ) 3 , d). Clearly, both polarization pro¬ 
files are strikingly different which in turn causes a loss 
of fidelity when the blockade radius is smaller than the 
system length. We verify this by numerically calculating 
the fidelity as a function of the blockade radius. The data 
is displayed in Fig. [ 4 ^ 1 , together with the corresponding 
transmission T. As anticipated, the fidelity decreases 
significantly below unity when R^ is decreased with re¬ 
spect to L. Note, that the transmission is close to zero 
throughout. 

A fidelity smaller than unity directly indicates the for¬ 
mation of a mixed state after the photon scattering. The 


0 0.1 0.2 0.3 0.4 


<7 

<D 

a 
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FIG. 3. (a,b) Squared modulus of the polarization \P(z,t)\ 2 
and (c,d) the time-integrated intensity I p (z) for R\> = L/2 
and g = IOOO 7 for two different positions of the gate atom. 
The gate atom (green circle) is located at Zi — 0 in panels 
(a) and (b) and at Zi — L in panels (c) and (d). The dashed 
line marks the blockade radius with respect to the gate atom 
position. 


final state density matrix is 

Pf~^Z A * e[ik(Zi ~ Zk)] \Zi)(Zk\. 

jk 


The final state can only be pure when \Aj k \ = 1 and 
hence F = 1. The formation of a mixed state is a conse¬ 
quence of the actual measurement of the gate atom po¬ 
sition [42] which is performed by the photon scattering: 
When R\) < L one in principle gains information on the 
position of the gate atom since the spatial uncertainty 
of its wave function is reduced from L to the blockade 
region. The final state is then a mixture of all states 
compatible with this additional information. 


V. ANALYTICAL RESULTS FOR COHERENT 
PHOTON SWITCHS 


In the remainder of the paper we will focus on the 
case of a coherent photon switch, i.e. R\> > L. Here the 
expression for the susceptibility of the medium simpli¬ 
fies to that of an ensemble of two-level atoms, Xj( z ) ~ 
— g/(uj + 27 / 2 ) which permits the derivation of analytical 
results. For a narrow band width pulse we can derive 
explicit solutions to Eq. 0 that have no dependence on 
the position of the gate atom m- For example, the po- 


















5 



FIG. 4. (a) Fidelity (circles) and transmission (squares) as 
a function of the blockade radius Rb with g = IOOO 7 . The 
transmission is already negligible when Rb/L > 0.4, but the 
fidelity approaches unity only when Rb ~ L. (b) Intensity 
I P (z) in the strong blockade regime (Rb = L) for g = IOOO 7 
(square) and g = 5127 (circle). The squares and circles are 
numerical data. The solid and dotted curves are the analytical 
results obtained from Eq. <©• 


larization P(z , t) is given by 


P(z,t) « ig\ 


/v^T 

cy 2 r 2 

- 47 (cT + z 0 ) 

4 g 2 z 

V 2c exp 


8c 

07 


xEc 


C 7 3 r 2 — 27 2 (cT + zq) — Sg 2 z 


2v / 2ct7 2 


(6) 


where Ec(x) is the complementary error function. The 
corresponding time-integrated profile I p (z) agrees per¬ 
fectly with the numerical result from Eq. Q (see Fig. Bo- 
The transmission T is given by 


T 



1 + Er 


L - z 0 

CT 



(7) 


where Er(x) is the error function and a = 4g 2 L/c~/ is 
the optical depth of a resonant two-level medium. The 
excellent agreement between the analytical and numeri¬ 
cal calculation is shown in Fig. It- Neglecting the finite 
band width of the photon pulse, i.e. when all the fre¬ 
quency components are in the absorption window, Eq. 0 
reduces to the well-known form T « e m- 

Finally, the fidelity can be expressed as a function of 
the optical depth and pulse band width 


F& (1 — e~ a ) 


2 12 

q/2'7—2 /y 4 /y -4 


( 8 ) 


This shows that indeed a small band width is a require¬ 
ment for reaching a large fidelity. For example, the trans¬ 
mission is negligible (T ^ 8 x 10 -4 ) when yr = 5 and 
g = IOOO 7 according to the data in Fig.[2^. However, the 
fidelity is below unity (F = 0.94) due to non-negligible 
contributions from the terms accounting for the finite 
band width. In the limit of very long pulses one finds 
F ~ 1 — T and thus the fidelity is solely determined by 
the transmission. 


VI. SUMMARY 

In summary, we have studied the coherence of a Ryd¬ 
berg spinwave in the operation of a signle photon switch. 
The current study is limited to a single gate atom and 
an incoming single-photon pulse, which permits the de¬ 
scription of multi-photon scattering, however, only if the 
photons enter the switch sequentially. Addressing this 
limitation and extending the discussion to correlated and 
entangled photon pulses that fall in the operation regime 
of single photon transistors will be subject to future stud¬ 
ies. 
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Appendix A: Details of the analytical calculation 


Here we will show how to obtain the analytical solution 
to the coupled equations (3) in the main text in the strong 
blockade regime (R^ > L) and for narrow band pulses. 

In the frequency domain, the solution to equations (3) 
is given by, 


£(z, id) = exp 


iuj 

z T 
c 


igz 


-x 


£o M, (Al) 


P{z,ui) = x£{z,ui), (A2) 

S(z,w) = -(A 3 ) 

UJ — V 


Due to the strong blockade condition, we have removed 
the dependence of £(z, cu), P(z,uj), and S(z,uj) on the 
gate atom index j and replaced the susceptibility by the 
one corresponding to two-level atoms, % = — g/(u+i'y/2). 
Moreover, we set S(z, uS) ~ 0 , which is a go od a pproxima- 
tion as Q/[u — V(z)\ ~ tl/V(z ) « 0 in Eq. ( |A3| ). Our aim 
is to obtain analytical expressions of £(z, t ) and P(z,t)- 
Applying the inverse Fourier transform on the both 
sides of Eqns. ( |A 1 | ) and (A 2 ), we obtain the formal solu¬ 
tion for £(z,t) and P(z, t) in time domain , 


/ oo 

du£ 0 (u)e- iujT+i ^ x , (A4) 

-OO 

/ OO 

dco X £o(u})e- iuT+i ^ x . (A5) 

-OO 
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The integration is in general difficult to carry out analyt- ity in powers of u;, 
ically due to the complicated form of the susceptibility. 

We overcome this difficulty by expanding the susceptibil- _ 


2 uj Aiuj 2 

- 2 + 

7 i 


(A 6 ) 


First let us calculate the approximate solution for £(z,t). To carry out analytical calculations and at the same 
time tak e into account contributions due to the finite band width, we will keep terms up to the second order of uj in 
Eq. (A 6 ). This yields the solution for £(z,t) 


£{z,t) = 


1 


■vV^crCO) 


exp 


2 r* 

cy 


1 


2 e(z)i 


t — 


7 2 - 4# 2 z 


fo 

c 


(A7) 


with £(z) = yFWW. For the current problem, we always have £(z) > 0 as the photon travelling time 
through the medium is the shortest time scale. For example, L/c'zz 3 x 10 “ 14 second for L = 10/xm. 

With the solution for £(z,t), we can calculate the transmission T = J 0 °° dt\£(L,t)\ 2 / J 0 °° dt\£ (0, £)| 2 . We need to 
carry out the respective two integrals over time at z = 0 and z = L. This can be done analytically, 

fOO 1 _ rJ Zp \ 

J q dt\£(0,t)\ 2 = -(A 8 ) 

and 


/ 


dt\£(L, t)\ 2 = 


2 c£(L) 


1 + Er 


/ L — zo a 
\ct£(L) 7 t£(L) 


2c£(L) 


1 + Er ( — 


zo 


a 


ct£(L) 7 t£(L) 


(A9) 


This leads to the analytical form of the transmission (7) in the main text. 

With the analytical solution for £(z,t) at hand, there are two ways to calculate P(z,t). We can directly calculate 
P(z,t) from Eq. (3b) by inserting the solution ( |A7| ) and S(z,t) = 0 . This yields the linear response of the medium to 
the photon electric field, 


P(z,t)=ig f dre l(yt t ^ 2 £(z,t) 

J —GO 

The integration over time can be carried out analytically, which gives 
P(z,t) 


W 


0 TT 


2 c 


exp 


cy 2 r 2 — 4q (cT + zq) 6g 2 z 


8 c 


07 


1 + Er 


3 £ 2 (z) cT + zo 

-7 T + 


4v/2£(^) y/2€(z)cr 


(A 10 ) 


(All) 


However it is difficult to calculate the fidelity from Eq. (All) due to the presence of the error function. We thus 
calculate P(z, t) a lternatively using the Fourier transform method. We note that the susceptibility x appears at two 
places in Eq. (A5): one in front of £(cz) and another one in the exponential function. In order to obtain an analytical 
result, we will expand the former susceptibility up to the second order of uj while the latter up to the linear order. 
After performing the inverse Fourier transform, we obtain the expression for P(z,t), 

P(M) = 32 ^ 4^9 {c 2 7 6 t 4 + 8 cz 7 3 tV - 4 c 2 7 3 r 2 + 2cy 5 T 2 (ct - z + z 0 ) + 4 [ 7 2 (ct - z + z 0 ) + 4g 2 z] 2 ) .(A12) 


Using Eq. (A 12 ) the fidelity can be calculated analyt¬ 
ically, 


(l-e"“) 


1 - 


y 2 r 2 


+ 


12 


7 4 r 4 


(A13) 
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